Superfluid Weight, Free Carrier Density, and Specific Heat 
of the d = 3 t J Model at Finite Temperatures 
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The superfluid weight, free carrier density, and specific heat of the three-dimensional tJ model 
are calculated by renormalization-group theory. We find that optimal hole doping for superfluidity 
occurs in the electron density range of (m) w 0.63-0.68, where the superfluid weight n s /m* reaches 
a local maximum. This density range is within the novel r phase, where the electron hopping 
strength renormalizes to infinity, the system remains partially filled at all length scales, and the 
electron-hopping expectation value remains distinctively non-zero at all length scales. The calculated 
superfluid weight drops off sharply in the overdoped region. Under hole doping, the calculated 
density of free carriers increases until optimal doping and remains approximately constant in the 
overdoped region, as seen experimentally in high-T c materials. Furthermore, from calculation of the 
specific heat coefficient 7, we see clear evidence of a gap in the excitation spectrum for the r phase. 
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I. INTRODUCTION 

The variation of the superfluid number density n s with 
temperature and carrier doping is of fundamental im- 
portance in describing the unique properties of the su- 
perconducting state in high-Tc cuprates. Experimen- 
tally, muon-spin-rotation techniques are used to deter- 
mine the closely related quantity n s /m* (also known as 
the superfluid weight), where m* is the effective mass 
of the carriers in the superfluid. In the undcrdoped re- 
gion of high-r c materials, n s /m* increases with doping, 
and the low-temperature superfluid weight is correlated 
with T c .0, 0] As the materials are doped past the op- 
timal value (where T, is the highest), n s /m* peaks and 
rapidly decreases. [!| |j, 0| The decrease in n s /m* is sur- 
prising since the total density of free carriers saturates 
at optimal doping and remains almost constant in the 
overdoped region. jfj] By contrast, in a conventional su- 
perconductor, described by BCS theory, these two quan- 
tities have the same doping dependence. 

The tJ model is a promising starting point in un- 
derstanding these properties of cuprate superconductors. 
Mean-field U(l) and SU (2) slave-boson theories of the t J 
Hamiltonian have reproduced some aspects of the dop- 
ing and temperature dependences of n s /m*.0,l3 More 
direct, unbiased numerical techniques applied to a 4 x 4 
tJ cluster have observed a large peak in n s /m* in the 
same region where pairing correlations indicate a super- 
conducting ground state. ;9| A general limitation of these 
types of studies is that no finite-cluster approach can 
unambiguously identify phase transitions in the system, 
or exhibit the non-analytic behavior of thermodynamic 
quantities at these transitions. 

Alternatively, the physics of the bulk model can 
be studied through the position-space renormalization- 
group method, which has been used to determine the 
phase structure and thermodynamic properties of the 
tJ and Hubbard models at finite temperatures. [Tol ITU 



IT2I Il3| In particular, Falicov and Berker's calculation 
for the t J model in d = 3 with the realistic coupling 
J/t = 0.444 produced a rich, multicritical phase dia- 
gram [Icllllj . with a novel low-temperature phase (called 
"t" ) for 30 — 40% hole doping where the electron hop- 
ping strength in the Hamiltonian renormalizes to infinity 
under repeated scale changes, while the system remains 
partially filled. This is the possible signature of a super- 
conducting phase, and it is notable that a similar phase 
was also observed in the d — 3 Hubbard model. [T2I IT^ | 

Our present study further develops this 
renormalization-group method, to yield the super- 
fluid weight of the t J model as a function of temperature 
and hole doping. Our approach reproduces phenomeno- 
logical features of high-T c materials. In particular we 
find that optimal doping is located in the vicinity of 
the t phase, where n s /m* peaks and then is sharply 
reduced with overdoping. Moreover, we also find that 
the density of free carriers increases until optimal 
doping, and saturates in overdoped region. These results 
suggest that the r phase might indeed correspond to the 
superconducting phase in cuprates. Further supporting 
this idea, we present specific heat calculations that show 
clear evidence of a gap in the quasiparticle spectrum for 
the r phase. 



II. THE tJ HAMILTONIAN 



We consider a d-dimensional hypercubic lattice where 
the t J model for electron conduction is defined by the 
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Hamiltonian 



H = P 



C i<T C j<T 



i 3 , (1) 



where and are creation and annihilation operators, 
obeying anticommutation rules, for an electron with spin 
a = | or J, at lattice site i, rii a = cj^c^, n, = + 

arc the number operators, and Si = J2acr> c t a s ^cr' c i^' * s 
the single-site spin operator, with s the vector of Pauli 
spin matrices. The entire Hamiltonian is sandwiched be- 
tween projection operators P — Yli(^ ~ n il n tf)> which 
project out states with doubly-occupied sites. The inter- 
action constants t, J, V describe the following physical 
features: electron hopping (t), a nearest- neighbor anti- 
ferromagnetic coupling (J > 0), and a nearest-neighbor 
interaction (V). The standard tJ Hamiltonian is a spe- 
cial case of Eq. with V/J = 1/4. For convenience, we 
introduce dimensionless interaction constants t, J, V, /z, 
and rearrange the jl chemical potential term to group 
the Hamiltonian into a single lattice summation: 

-p h = £ p [-* £ (4<> + sW) 



— JSi ■ Sj + VniTij + (x(rii + nj) 



P (2) 



Here j3 = 1/fcsT, so that the interaction constants are 
related by t = (3t, J = (3 J , V = (3V , n = (3p,/2d. The 
total Hamiltonian is now written as a sum of pair Hamil- 
tonians —/3H(i,j). The sum over nearest-neighbor sites 
is taken so that the position of site j is rj = Vi + ajt, 
where is one of the d lattice vectors. Since changing 
the sign of t is equivalent to redefining the phase at every 
other site in the system, we shall choose t > with no 
loss of generality. The effective temperature variable will 
be 1/t = ksT, where we have taken t = 1 as the unit of 
energy. 

In order to study the superfluid weight, we introduce 
periodic boundary conditions, by considering the system 
as a ring in each axis direction threaded by a magnetic 
flux. We choose the vector potential A associated with 
the flux to have components Aj\fd along each axis, so 
that the pair Hamiltonian becomes 



f3H(i,j)=P 



-*£ 



6 "^Cj^Cjcr -f" C ^Cj^Ci, 



JSi • Sj + VriiUj + [i[ni + rtj 



P, (3) 



where <fi — aA/y/d and a is the lattice spacing. For sim- 
plicity, we have adopted units so that ft = c = e = 1. In 
the presence of the applied phase twist <j>, the superfluid 
weight n s /m* is related to the curvature of the total free 
energy F near 



n s 1 
= ^ hm 



d 2 F 



No 2 a^o dA 2 



1 d 2 F 



(4) 



where N — ► oo is the total number of lattice sites. In 
Sec. HIE we shall show how this quantity can be calcu- 
lated from the renormalization-group transformation de- 
veloped below. 



III. RENORMALIZATION-GROUP 
TRANSFORMATION 

A. Recursion Relations 

The position-space renormalization-group method 
used here starts with an approximate decimation in 
d = 1, which is then generalized to higher dimensions 
by the Migdal-Kadanoff procedure [HIE]]. In d = 1, the 
Hamiltonian of Eq. J5J takes the form: 



-(3H = Y,{-(3H(i,i + l)} 



(5) 



where i = 1,2, 3, . . .. The decimation consists of finding a 
thermodynamically equivalent system, described by the 
Hamiltonian —f3'H', which depends only on the states of 
the odd-numbered sites. Since the quantum operators in 
the Hamiltonian do not commute, an exact decimation 
even in one dimension is not possible. We can carry out 
an approximate decimation as follows |l6lll7|: 

Tr p-P H -Tr P J2i{-0H{i,i+l)} 

J- 1 even J- 1 nvnn c- 



1 even ^ 
even 



(6) 



TT e -p'H'(i-i,i+i) 

i 



Here —fl'H 1 is the Hamiltonian for the renormalized 
system, and Tr eV en is a trace over the degrees of free- 
dom at all even-numbered sites. In the two approxi- 
mate steps, marked by ~ in Eq. ©, we ignore the non- 
commutation of operators separated beyond three con- 
secutive sites of the unrenormalized system (conversely, 
this means that anticommutation rules are taken into 
account within three consecutive sites at all successive 
length scales, as the renormalization-group procedure is 
repeated). These two steps involve the same approxima- 
tion but in opposite directions, which gives some mutual 



3 



compensation. Earlier studies of quantum spin systems 
have shown the success of this a ppr oximation at predict- 
ing finite-temperature behavior. Il7j 

The renormalization-group mapping can be extracted 
from the third and fourth lines of Eq.©: 



,-0'H'{i,k) _ 



Trje 



-0H(i,j)-0H(J,k) 



(7) 



n 


V 


s 


m 3 


Two-site basis states 





+ 








\<f>i) = | o o) 


1 


+ 


1/2 


1/2 


|<M = ;75{|t°> + l°m 


1 




1/2 


1/2 


!<M = ^{|t°>-|oT)} 


2 










I0 6 ) = -^{| U> - 1 XT)} 


2 


+ 


1 


1 


\<h) = | TT) 


2 


+ 


1 





109) = ^{| U) + I IT» 



where k are three consecutive sites of the unrenormal- 
ized system. The operators —fi'H'{i, k) and —f3H(i,j) — 
(3H(j, k) act on the space of two-site and three-site states 
respectively, so that, in terms of matrix elements, 



[UiV k \e 



-0'H'(i,k)\ 



H u kl 



(8) 



where Ui 7 Wj,Vk,Ui,v k are single-site state variables. 
Eq.(JH| is the contraction of a 27 x 27 matrix on the right 
into a 9 x 9 matrix on the left. We block-diagonalize 
the left and right sides of Eq. (JSJl by choosing basis states 
which are the eigenstates of total particle number, total 
spin magnitude, total spin z-component, and parity. We 
denote the set of 9 two-site eigenstates by {|0 P )} and the 
set of 27 three-site eigenstates by {l^?)}: an d list them 
in Tables I and II. Eq. (JHJ) is rewritten as 



-0'H'(i,k) 



\<h) 



J2 E<^i u * ,, *)< u * ti, j ,, *i^)^i e " w<J) " /3 ^ (i ' fc) i^) 



u,v,u, q,q 
v.w 



{tl}il\uiWjV k ){uiV k \(f)p) . (9) 

Eq. J5J| yields six independent elements for the matrix 
\ e -0 h ^ which we label 7 p as follows: 



7 P 
7o 



-p'H'(i,k) 



\<t> p ) for p = 1,2, 4, 6, 7, 



He -0'H' {i ,k)^ A) 



(10) 



To calculate the 7 P , we determine the matrix elements 
of —(3H(i,j) — /3H(j,k) in the three-site basis {tp q }, as 
listed in Table III, and exponentiate the matrix blocks to 
find the elements (\j} q \e~ l3H ^^~ l3H ^^\'il' q ) which enter 
on the right-hand side of Eq. |JSJ| ■ In this way the 7 P are 
functions of the interaction constants in the unrenormal- 
ized Hamiltonian, 7 P = 7 P (t, 0, J, V, /i). 

Since the matrix (0 p |e - ^ H ^ % ^\4>p) is determined by 
six independent elements 7 P , the renormalized pair 
Hamiltonian —@'H'(i,k) involves six interaction con- 
stants, namely those of the original types of interactions 
and an additive constant: 



(3'H'(i,k) = P 



e ^ c i(J Cj a + e 



J' Si ■ Sj + Vrnnj + n'(m +Uj) + G' 



P, (11) 



TABLE I: The two-site basis states, with the corresponding 
particle number (n), parity (p), total spin (s), and total spin 
z-component (m s ) quantum numbers. The states \<f>3), |0s), 
and |08), are obtained by spin reversal from |</>2), 104), and 
|07), respectively. 



n 


P 


s 


m a 


Three-site basis states 





+ 








k/'i) = | o oo) 


1 


+ 


1/2 


1/2 


|Va> = l°t°>, hfe) = ^{|T°°) + |° °T)} 


1 




1/2 


1/2 


kM = ^{|T°°)-|°°T>} 


2 


+ 








hM = 3{|tl °)-UT °)-|°tl> + MT)} 


2 










liM = |{| U °>-UT °) + l°TI)-|°IT)}, 
|^io) = ^{|t°l)-U°t)} 


2 


+ 


1 


1 


hM = [t°t>. hM = -^{| TT °) + |oTt)} 


2 


+ 


1 





\*h*) = i{\U °) + |IT °) + l°T4) + l°IT)}, 
hM = 73{|W> + |W>} 


2 




1 


1 


liM = -75 TT TT °) - 1° TT)} 


2 




1 





|-0is) = HI TT °) + U?°)-|°U)-|°TT)} 


3 


+ 


1/2 


1/2 


IVao) = ^{2| TTT) - 1 TTT) - 1 TtT)} 


3 




1/2 


1/2 


IV22) = ^{| TU) - 1 m» 


3 


+ 


3/2 


3/2 


1^24) = | TTT) 


3 


+ 


3/2 


1/2 


1^25) = -75 iT TIT) + I TTT) + 1 TTT)} 



TABLE II: The three-site basis states, with the correspond- 
ing particle number (n), parity (p), total spin (s), and total 
spin z-component (m a ) quantum numbers. The states |04_s), 
107), 1015-ie), |0ig), |02i), |023), 1026-27), are obtained by 
spin reversal from |02-3), 106), </>n-i2> , |0rr), |02o), 1022), 
1024-25), respectively. 



The matrix elements of —[3'H'(i,k) in the {0 P } ba- 
sis are shown in Table IV. Exponentiating this matrix, 
we can solve for the renormalized interaction constants 
(f, 0', J', V', fJ,', G') in terms of the 7 P : 

t' = sign (74 - 72) coslT 1 ( j 2 ^ 4 ^) > 



b' = tan 



2 Im.70 
74 - 72 



V = 1 {ln(7i 4 767 7 3 ) - Zv} , 



J' = In—, 

76 

jjf = v- In 71, 



G' = ln 7l , 



(12) 



where 



-In (7274 - 7o7o) 



The approximate d = 1 decimation contained in 
Eqs. (|9"|)- l|12[) can be expressed as a mapping of a Hamil- 
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Vi 


01 










■02 


■03 


06 






— \Tlt COS 


2£ sin 




— \/2t cos 




iAi 


^6 


— i \/2t sin 


-iAi 









^8 


•09 


•010 


•08 


+ y + 3m 


-iA 3 


— iy/2t sin 


•09 




-f J+ V + 


— cos 


V>10 


i^/2t sin 


— \/2t cos 


2p 





0>11 


</>12 


•017 


•011 


2m 


— cos 


z v^i sin 


•012 


— COS 


ij + y + 3 M 


iA 2 


•017 


— iy/2t sin 


-iA 2 


i ./ + V + 3m 





•013 


•014 


•018 


•013 


i,7 + V + 3^ 


— v^i cos 


iA 4 


•014 


— cos 


2jU 


sin 


0>18 


-1A4 


— z'V2t sin 


i .7 + V + 3/j 





•020 


"020 


-.7 + 2V + 4m 





■022 


•022 


2V + 4m 





0124 


0124 







■025 


>/>25 


^./ + 2V + 4m 



TABLE III: Diagonal matrix blocks of the unrenormalized 
three-site Hamiltonian —/3H(i,j) — f3H(j,k). The Hamilto- 
nian being invariant under spin-reversal, the spin-flipped ma- 
trix elements are not shown. The additive constant contribu- 
tion 2G, occurring at the diagonal terms, is also not shown. 
The additional A, terms, which are not part of the original 
three-site Hamiltonian, are explained in Sec.IIIC. 



decimation described above to the middle site between 
the two consecutive segments, giving the renormalized 
nearest-neighbor couplings for the points in the new lat- 
tice. We compensate for the interactions that are ignored 
in the original lattice by multiplying by a factor of 6 d_1 
the interactions after the decimation, 6 = 2 being the 
length rescaling factor. Thus, the renormalization-group 
transformation of Eq. l|13fl generalizes, for d > 1, to 



K' = 6 d - 1 R(K) 



(14) 



B. 



Renormalization-Group Transformation in the 
Presence of Magnetic Flux 



In order to correctly model the response of the sys- 
tem to an applied magnetic flux, the renormalization- 
group approximation described in the last two sections 
needs to be extended. To see this, we first review the 
formalism for calculating thermodynamic densities from 
the renormalization-group flows. [20| Conjugate to each 
interaction K a of K = {K a }, there is a density M a (e.g., 
kinetic energy, electron density), 



1 dlnZ 
~Nd 8K a 



(15) 



where Z(K) is the partition function. We can re- 
late the densities at the two consecutive points along a 
renormalization-group trajectory by 





01 


02 


04 


06 


07 


09 


01 


G' 
















— t' COS 0' + 


it' sin 0' 









02 




m' + G' 










—it sin 


£' cos 0' + 








04 




M' + G' 








06 








-%J'+Vi + 
2m' + G' 






07 











2m' + G' 




09 












2m' + G' 



TABLE IV: Block-diagonal matrix of the renormalized two- 
site Hamiltonian —f3'H'(i,k). The Hamiltonian being invari- 
ant under spin-reversal, the spin-flipped matrix elements are 
not shown. 



tonian with interaction constants K = {G,t, J,V, fi,6} 
onto another Hamiltonian with interactions constants 

K' = R(K) . (13) 

The Migdal-Kadanoff procedure [TsL ITflj is used to con- 
struct the renormalization-group transformation for d > 
1. We ignore a subset of the nearest-neighbor interactions 
in the d-dimensional hypercubic lattice, leaving behind a 
new d-dimensional hypercubic lattice where each point 
is connected to its neighbor by two consecutive nearest- 
neighbor segments of the original lattice. We apply the 



b- d M' p Tp 



where 



T, 



(3a 



9K' p 
8K a 



(16) 



with summation over repeated indices implied. At a fixed 
point of the renormalization-group transformation, corre- 
sponding to a phase transition or a phase sink, the densi- 
ties M a = M' a = M* are the left eigenvector with eigen- 
value b d of the recursion matrix T evaluated at the fixed 
point. The densities at the starting point of the trajec- 
tory (the actual physical system) are computed by iter- 
ating Eq. i|16|) until a fixed point is effectively reached. If 
T( fc ) is the recursion matrix of the (fc)th renormalization- 
group iteration, then for large k, we can express the den- 
sities of the actual system M as 



M ~ b- kd M* ■ [T (fc) ] • [T^ 1 '] [T 



Ml 



(17) 



The renormalization-group transformation incorpo- 
rated in Eqs. ©"031 g ives 



(18) 



for all 6. The 6x6 recursion matrix T will then have 



df 


83' 


8V 


= ^ = 


86' 
~8J 


~ 90 ~ 
_ d $ _ 

" w ~ 


~8§ 
86' 
8\i 


86 ' 
86' 

— n — — — 2 

_U ' 86 < 



5 



the form 



I u at 



dt 
dt 



u at 



\ 



ac 
a^L 
at' 

dfj. 



a^ 



(19) 



at every step in the flow. This leads to 

MS =0 and — In Z = M 6 = , (20) 
ocp 

for all points of the phase diagram. This superfluid 
weight of zero for all temperatures and electronic den- 
sities is clearly due to the oversimplification in our initial 
approximation. 

The source of the problem is the three-site cluster ap- 
proximation used in deriving the recursion relations. In 
modifying the original approximation scheme, we seek to 
incorporate the effect of the non-commutations extending 
beyond the three-site cluster. Turning to the matrix ele- 
ments of —/3H(i, j) — (3H(j, k) listed in Table III, we note 
the terms A,*, i = 1, ... ,4. Using the original Hamilto- 
nian of Eq. (|3J restricted to the three-cluster, the matrix 
elements involving these terms are all zero. However, 
non-commutativity extending beyond the three-cluster 
makes, as we see below, these matrix elements non-zero. 

We can estimate the magnitude of the matrix elements 
Aj by considering a five-site cluster, described by Hamil- 
tonian -0H(h, i) - f3H(i, j) - /3H(j, k) - (3H(k, I), where 
(h, i, j, k, I) are consecutive sites. In the spirit of Eq. iJHJ, 
we generate effective couplings for the three-cluster by 
tracing over the degrees of freedom at the outside sites 
in the five-cluster, 

(m Vj w k \ e -PH(i,j,k) i - . ^ = 

■ \thUiVjWkXi) , (21) 



where the subscripted variables refer to single-site states. 
From the above equation, we can extract the ma- 
trix elements of an effective three-cluster Hamiltonian 
-j3H(i,j,k). Eq. (J2TJ is the contraction of a 243 x 243 
matrix on the right-hand side into a 27 x 27 matrix on 
the left. We simplify our task by using the {ip p } basis 
on the left, and choosing an appropriate five-site basis to 
block-diagonalize the 243 x 243 right-hand matrix. 

Since —/3H(i,j,k) is derived from the decimation of 
a five-cluster, it will have a more general form than 
—j3H{i,j) — (3H(j,k), and approximately reflect the ef- 
fect of the three-cluster non-commutations with the ex- 
ternal sites. However our approximation scheme must 
also satisfy an important constraint: the ef> — > limit 
should yield the same renormalization-group transforma- 
tion used in earlier studies of the tJ model [Tfl Hl| . To 



achieve this, we modify only a subset of the matrix ele- 
ments of —(3H(i,j) — /3H(j,k), namely those which are 
zero in the original scheme when <f> ^ 0, but whose cor- 
responding elements in —(3H(i,j,k) are non-zero: 



Ai = sign(cj)t)\{tp 3 \0H(i, j, k)\tp 6 }\, 
A 2 =sign(#)|(V>i 2 |^(z,j,fc)|Vi7)|, 
A 3 =sign(0t)|^ 8 |^(i,i ! fc)|V9)|, 
A 4 = sign(#) | (*p 13 0H(i, j, k)\ip 18 ) | . 



(22) 



The sign(4>t) prefactors guarantee that couplings between 
the same types of three-cluster states have the same 
sign. For example, \tp2) and ^3) share the same n, p, 
s, and m s quantum numbers, as can be seen from Ta- 
ble II. A nonzero <j) couples \ifj 2 ) to \ipe), a state with 
the same n, s, and m s , but opposite parity. From the 
second block in Table III, the associated matrix ele- 
ment is (ip2 \ ■ • • |"06) = iv2i sin(<^) . The Ai elements in 
that block have an analogous role, coupling \ip 3 ) to j^)- 
The prefactor in the Ai expression of Eq. (|22|l sets the 
sign of the element (?/> 3 | • • • |V>6) = *Ai to equal that of 
(-02 1 • • • I'ipe)- Since our calculations are all done for small 
<f>, sign(sin0) = sign(0). Similar reasoning applies to the 
prefactors of the other A^ elements. 

Through Eq. (|21|1 . the A^ are functions of the in- 
teractions strengths in the unrenormalized Hamiltonian, 
Ai — Ai(t, J, V, /1, 4>). They scale like 4> for small </>, and 
duly vanish in the limit <f> — > 0. As will be explained 
in Sec. HIE, finding the superfluid weight involves calcu- 
lating a thermodynamic density in the cf> — > limit, so 
we shall be working in the regime where the A^ are van- 
ishingly small. The result of the extended calculation, 
taking into account the quantum mechanical backflow 
into the three-cluster, is that Eqs. Ijl8|l no longer hold, 
dhiZ/d(j) ^ in general, and we obtain interesting non- 
trivial results for n s /m*. 



C. Calculation of the Superfluid Weight 

The superfluid weight of Eq. (@J is expressed as a 
derivative of the total free energy F — F(n, T, 0), where 
n = (rii) is the electron density. In terms of the conjugate 
current 



j(n,T, 
Eq. becomes 



1 dF 
~Nd 



n,T 



(n,T) = lim ■— ■ 

m* 0^0 dd) 



n,T 



In terms of the grand potential f2(At, T, 4>) 



,( M ,T^) = — — 



(23) 
(24) 

(25) 



6 



and 



n(jjt, T, 4>) 



2Nd ~dn~ 



(26) 



Relating the partial derivatives of j with respect to 
through 





dj _ dj 
d<f) T dn 


dn 

0,T d( i> 


H,T 


dj 
H - 

d<P 


n,T 


(27) 


and using 


the Maxwell relation 


dn 
d<f> 




/3 dj 
2 dn 


4>;T 


dj 


= _ii_dj_ 


dj_ 




dj 

dcj) 




(28) 


dcj) 


rp 2 dn 




<t>,T 


n,T 



The current j is zero when <fi = 0, so that the first term on 
the right-hand side above is also zero in the limit <p — > 0, 



and we find that lim0_>o 



Thus 



Eq. 124|) can be equivalently written as 



— ^-(/i,T) = lim — 
m* ^->o dm 



u,T 



Nd 0™ d<h 2 



ft,T 



1 ,. d 2 \nZ 
hm 



(3Nd <l>^o 



(29) 



ti,T 



This is the form we shall use when calculating the super- 
fluid weights. 



u 
ft 
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FIG. 1: Phase diagram for the d — 3 tj model with J/t = 
0.444, 4> — 0, in temperature versus electron density, jl^ The 
antiferromagnetic (A), dense disordered (D), dilute disordered 
(d), and r phases are seen. The second-order phase bound- 
aries are drawn with full curves. The coexistence boundaries 
of first-order transitions are drawn with dotted curves, with 
the unmarked areas inside corresponding to coexistence re- 
gions of the two phases at either side. The dashed lines are 
not phase transitions, but disorder lines between the dilute 
disordered and dense disordered phases. 
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IV. RESULTS 



TABLE V: Expectation values at the phase-sink fixed points. 



A. Global Phase Diagram for d — 3 

Each sink, or completely stable fixed point of the 
rcnormalization-group flows, corresponds to a thermo- 
dynamic phase, and we find the global phase diagram by 
determining the basin of attraction for every sink |21| . 
Flows that start at the boundaries between phases have 
their own fixed points, distinguished from phase sinks by 
having at least one unstable direction. Analysis of these 
fixed points determines whether the phase transition is 
first- or second-order. As explained in Sec.lllC, the ther- 
modynamic densities, which are the expectation values 
of operators occurring in the Hamiltonian, can also be 
calculated from the renormalization- group flows. In par- 
ticular, we determine the single-site electron density (n*). 
For the coupling J/t = 0.444 and <j> = 0, the phase di- 
agram in terms of (rij) and temperature l/t is shown in 

Fig. i nana. 

The nature of the various phases is epitomized by the 
thermodynamic densities M* calculated at each phase 
sink (Table V), which underpin the calculation of den- 
sities throughout their respective phases (Eq. I|17|)l. We 
summarize the phase properties below (for a more de- 
tailed discussion, see \wi fflf h 



Dilute disordered phase (d): The electron density 
(ni) = at the sink and, as a result, the (n,) calculated 
inside this phase are low. 

Dense disordered phase (D): The electron density 
(rii) = I at the sink and, as a result, the (nA calculated 
inside this phase are close to 1. 

Antiferromagnetic phase (A): The electron den- 
sity (n^ = 1 at the sink, so that this phase is also 
densely filled. The nearest-neighbor spin-spin correla- 
tion (Si • Sj) = 1/4 at the sink. Two spins that are 
nearest neighbors at the sink are distant members of the 
same sublattice in the original cubic lattice. The non- 
zero value of the correlation function at the sink leads 
to (Si • Sj) < for nearest-neighbor sites of the original, 
unrenormalized system. 

r phase: This is a novel phase, characterized by 
partial-filling at the phase sink, (n,) = 2/3. It is the 
only phase where the electron hopping strength t does 
not renormalize to zero after repeated rescalings; instead, 
t — > oo at the sink. As a result, the expectation value 
of the electron hopping operator at the sink is non-zero, 

J2a{ c icr c 3<r + c ]<j c i<?) = -2/3. This property makes it a 
possible t J model analogue to the superconducting phase 
in high-T c materials. The superfluid weight and thermo- 
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dynamic results discussed below certainly support this 
idea. 

In the limit (rii) = 1, the system exhibits antiferro- 
magnetic order at low temperatures, as expected from 
the spin-spin coupling in the Hamiltonian. Upon hole 
doping, there is a competition between the A and D 
phases, which respectively minimize antiferromagnetic 
potential energy and hole kinetic energy. Note the ex- 
tent of the A phase near (rii) = 1, which persists only 
up to a small amount of hole doping 5 = 1 — (rii) < 
0.05. This feature is directly reminiscent of the an- 
tiferromagnetic phase in certain high-T c materials, for 
example L^-^Sr^CuC^ |22j . At intermediate dopings 
6 w 0.32 — 0.37, we have a low-temperature r phase, sur- 
rounded by islands of antifcrromagnetism. (When the 
hopping strength t increases under rescaling, this also 
lowers the free energy of antiferromagnetically long-range 
ordered states, which leads to these islands of A in the 
vicinity of the r phase. [T0|) At hole dopings S > 0.37, 
there is a transition to a dilute disordered phase, with a 
narrow region of first-order coexistence at lower temper- 
atures. 



B. Superfiuid Weight and Kinetic Energy 

Using the method of calculating thermodynamic densi- 
ties described in Sec. HIE, we determine (l/Nd)dlnZ/dcj> 
at small (j). Taking the numerical derivative of this quan- 
tity at (j> = 0, we find n s /m* through Eq. The su- 
perfiuid weight is plotted as a function of electron density 
in Fig. 2, along four different constant temperature cross- 
sections of the phase diagram. For comparison, we also 
show in the same figure the calculated average kinetic en- 
ergy per bond (K), where K = - J2 a (4a c j<r + c )a c i°) ■ 
(K) and the total weight of a\{uj,T), the real part of the 
optical conductivity, are related by the sum rule [2j| , 



J dua 1 (u;,T) = ^-(K) 



(30) 



In comparing the properties of the t J model to those of 
high-T c materials, we keep in mind that the t J Hamilto- 
nian describes a one-band model, so cannot account for 
interband transitions. For real materials, the full con- 
ductivity sum rule has the form 



i d,U>(Tl(LJ,T) 

Jo 



ire 2 n 
2m 



(31) 



where n is the total density of electrons and m is the free 
electron mass. The right-hand side of Eq. (|31(l is inde- 
pendent of electron-electron interactions, in contrast to 
the right-hand side of Eq. (j30|l . where (K) will vary with 
the interaction strengths in the Hamiltonian. The optical 
conductivity of actual materials incorporates both tran- 
sitions within the conduction band and those to higher 
bands, while the tJ model contains only the conduc- 
tion band. We can look at Eq. (|30|l as a partial sum 



rule |23j,|2J], which reflects the spectral weight of the free 
carriers in the conduction band. 

The experimental quantity we are interested in model- 
ing is the effective density of free carriers, which in actual 
materials is calculated from the low-frequency spectral 
weight psl ]. 



2m 



duj <ji (u>, T) 



(32) 



For high T c materials, the cutoff frequency is typically 
chosen around Tiloq ss 1 eV so as to include only intra- 
band transitions. For comparison with the t J model, we 
identify the right-hand side of Eq. I|30|l with 



7re 



(K) = 



ire 2 n flcc (T) 
2m 



The superfiuid weight satisfies the inequality 26 



<(K) 



(33) 



(34) 



which is obeyed in our results in Fig. 2. 

The superfiuid weight graphs at the sampled temper- 
atures show a clear bipartite structure, with a peak at 
low (rii), and another peak at high (rii) (which develops 
into two closely spaced peaks at lower temperatures). In 
between these is a region of low superfiuid weight, with 
a minimum near (rii) — 0.385, approximately indepen- 
dent of temperature. Looking at the nearest-neighbor 
density-density correlation (riirij) as shown in Fig. 3, we 
see that (rii) — 0.385 is also the electron density separat- 
ing two different regimes of the system: an extremely di- 
lute regime, where (riiiij) ~ 0, and a partially-to-densely 
filled regime, where (riirij) > 0. It is therefore useful to 
discuss the superfiuid weight and kinetic energy results 
in terms of these two regimes. 



1. Extremely dilute regime, (rii) ^ 0.385 

The system in this regime is a dilute gas of electrons. 
For low (rii), the kinetic energy per bond (K) ~ 2(rii), 
which follows if the density of free carriers is just the den- 
sity of electrons, nf rco = (rii), and the mass of the carriers 
m = 1/2. The interaction terms in the tJ Hamiltonian 
create an attractive potential of strength —J between 
electrons in singlet-states on neighboring sites. For a 
coupling J/t = 0.444, this attraction is too weak to form 
two-body bound states, but since we are in three dimen- 
sions, even a weak attractive potential is sufficient for 
the formation of an electron superfiuid at low tempera- 
tures H2i|28|. In fact, we see a peak in n s /m* develop 
around {rii) ~ 0.3-0.35, and this peak grows as the tem- 
perature is lowered from 1/t = 0.315 to 0.1. For low 
(rii), the superfiuid weight increases with electron den- 
sity and (K). The location of the peak in n s /rn* is just 
before (K) comes to its maximum and levels off. As the 



8 




0.2 0.4 0.6 0.8 
Electron density <«,-> 

FIG. 2: The superfluid weight n a /m* (solid line) and free 
carrier density n/ ree /m (dotted line) as a function of elec- 
tron density at four different values of temperature 1/t. The 
corresponding phases are indicated above the plots, and the 
location of phase boundaries marked by thin vertical lines. 
The symbol p refers to a region of forbidden densities due 
to the discontinuity at a first-order transition. The symbol 
L refers to a "lamellar" region where narrow slivers of the A 
and D phases alternate. 



density of free carriers saturates near (n^) ~ 0.385, there 
is a sharp drop in n s /m*, and (niUj) begins to increase 
from zero. At this density the physical characteristics of 
the system abruptly change, without however inducing a 
phase transition. 




0.2 0.4 0.6 0.8 
Electron density <«,-> 



1.0 



FIG. 3: The nearest-neighbor density-density correlation 
{rurij) as a function of electron density at two different values 
of temperature 1/t. The corresponding phases are indicated 
above the plots, and the location of phase boundaries marked 
by thin vertical lines. The symbol p refers to a region of 
forbidden densities due to the discontinuity at a first-order 
transition. The symbol L refers to a "lamellar" region where 
narrow slivers of the A and D phases alternate. 



2. Pariially-to- densely filled regime, (rij) > 0.385 

For intermediate densities (rii) i=s 0.385-0.63, the ki- 
netic energy (K) remains approximately constant. Near 
(iii) ~ 0.63, there is a phase transition to a densely filled 
phase (either D or ^4). We go from a physical picture 
where the carriers are electrons in a mostly empty back- 
ground to one where the carriers are holes moving in a 
mostly filled background. These holes condense into a 
superfluid at lower temperatures, and the peak in n s /m* 
occurs in the vicinity of the dilute-dense narrow first- 
order phase transition. For 1/t < 0.16, the maximum 
superfluid weight is reached inside the r phase. In the 
densely filled regime, (rii) > 0.63, the kinetic energy goes 
linearly as (K) ~ 2(1 — (rii)) = 2<5, as expected if the free 
carriers are holes. 

For hole-doped high-T c materials, the density of free 
carriers increases with 5 until the doping level optimal 
for superconductivity is reached, and remains approxi- 
mately constant in the overdoped regime. [6j The super- 
fluid weight, in contrast, peaks near optimal doping and 
sharply decreases with overdoping. These trends are re- 
produced in our numerical results, identifying, from our 
calculated n s /m* maxima, the optimal doping for the t J 
model as 5 ~ 0.32-0.37, the range of densities where the r 



9 



0.45 



0.3 - 



0.15 - 



0.3 - 



0.15 - 



Ja) 

***» \ 
* \ 

\\ 
v\ 

^s^^ 


1 

0.980 

— 0.850 - 

VI. 1 JU 

0.671 


(b) 


1* 




0.6605 




0.6550 - 




U.olMJU 


•••* 





0.1 0.2 0.3 

Temperature l/t 



0.4 



0.9995 
0.980 . 
0.950 
0.750 
0.700 ■ 
0.670 




0.2 0.3 0.4 0.5 
Temperature l/t 



0.7 



FIG. 4: The superfluid weight n s /m* as a function of tem- 
perature l/t for various electron densities (jii) indicated in 
the legends. Fig. 4(a) shows results in the range of small to 
optimal hole doping, while Fig. 4(b) shows results for hole 
overdoped systems. 



FIG. 5: The specific heat per bond C, in units of fa, as a 
function of temperature l/t for various electron densities (rii) 
(indicated in the legends). Fig. 5(a) shows results in the 
range of small to optimal hole doping, while Fig. 5(b) shows 
results for hole overdoped systems. The small discontinuities 
in the plot for {rii) — 0.655 reflect temperature ranges where 
that particular density does not appear because of the narrow 
first-order phase transition. 



phase occurs. Note that optimal doping for high-T c ma- 
terials is lower than this, typically around S = 0.15, and 
the closely spaced double-peak structure of n s /m* at low 
temperatures near optimal doping is not observed. On 
the other hand, our approximation for the d = itJ model 
is closer to experiment in this respect than earlier nu- 
merical studies of the t J model, which focused mostl y on 
finite-cluster techniques applied to the d = 2 system [2!j . 
In these earlier studies optimal doping is identified near 
{rii) = 0.5 on the basis of d-wave pairing correlations and 
the peak in the superfluid weight 9] . Also in these earlier 
studies, the kinetic energy has a maximum at {rii) = 0.5, 
but, unlike experiments, does not saturate with overdop- 

To complete the description of the superfluid weight 
in this regime, in Fig. 4 we show n s /m* as a func- 
tion of temperature l/t at various electron densities (rii). 
For systems with small to optimal hole dopings, shown 
in Fig. 4(a), there is a clear onset temperature near 
l/t ~ 0.2 below which the superfluid weight rises rapidly, 
until it levels off near zero temperature. This behav- 
ior is in good comparison with experimental results with 
YBa 2 Cu 3 6+a; |21|. As we move past optimal doping 
to the overdoped systems of Fig. 4(b), we see a marked 
change in behavior, with the low temperature n s jm* sup- 
pressed. 



C. Specific Heat 

Since the superfluid weight peaks inside the r phase 
at low temperatures, it is interesting to check whether 
the t region has any other general characteristics of a 
superconducting phase. We have added a magnetic field 
spin coupling term to the t J Hamiltonian and have shown 
that the r phase continues to exist when H ^ 0, up to 
a critical field H C (T), which decreases with increasing 
temperature and goes to zero at the temperature of the 
t phase boundary. In our present study, we look at the 
spectrum of excitations of the system through the specific 
heat per bond 



C(n,T) 



d(H(i,j)) 



dT 



(35) 



calculated for = 0. If the r phase corresponds to the 
superconducting phase in real materials, we should see 
evidence of a gap in the excitation spectrum. 

The results for C as a function of temperature 1 jt are 
plotted in Fig. 5 for a series of different electronic den- 
sities (rii). Starting at (rii) — 0.9995, the smallest hole 
doping shown in Fig. 5(a), we observe a broad peak 
around l/t ~ 0.33, corresponding to fc^T ~ 0.75 J. We 
can identify this peak with the thermal excitation of the 



10 



N 0.2 



d p, L-+- 


T 

1 


D 


A 


D 


A 




V 











1 1 ' ^ " — 1 ^ — 1 

0.5 0.6 0.7 0.8 0.9 1.0 



Electron density <n ( > 

FIG. 6: The specific heat coefficient 7 = C'/T, in units of 
k%, as a function of electron density (rn) at temperature 
1/t = 0.015. The corresponding phases are indicated above 
the plot, and the location of phase boundaries marked by 
thin vertical lines. The symbol p refers to a region of forbid- 
den densities due to the narrow discontinuity at a first-order 
transition. The symbol L refers to a "lamellar" region where 
narrow slivers of the A and D phases alternate. 

spin degrees of freedom. As we dope the system with 
holes, the weight under the curve at lower temperatures 
increases due to excitation of charge degrees of freedom. 
As we approach optimal doping, a second peak develops 
around 1/t ~ 0.2. Note that this approximately coincides 
with the onset temperature below which we see a dra- 
matic increase in n s /m* in Fig. 4. The spin-excitation 
peak is also enhanced for (rii) w 0.65-0.75, which is re- 
lated to the appearance of an antiferromagnetic island 
around 1/t ~ 0.3 in that density range. 

The peak at 1/t ~ 0.2 grows rapidly near optimal dop- 
ing, reminiscent of the specific heat anomaly of high- 
Tc materials [32I H3- For (m) = 0.655 we see the ap- 
pearance of two subsidiary peaks below the main one at 
1/t ~ 0.2. These smaller peaks may be related to the 
complicated lamellar structure of A and D regions above 
the r phase boundary. For temperatures 1/t < 0.16, in- 
side the t phase, the specific heat is strongly suppressed, 
reflecting the opening up of a gap in the excitation spec- 
trum. We can see this gap more directly by looking at 
the low-temperature limit of the specific heat. Quasipar- 
ticle excitations contribute a linear term to the specific 



heat C ~ 7T for small T. In Fig. 6, we plot 7 = C/T 
as a function of electron density at a low temperature, 
1/t = 0.015. The specific heat coefficient 7 ~ in the 
A phase near half-filling, but then grows with increasing 
hole doping. At the onset of the r phase a gap opens 
in the quasiparticle spectrum, 7 falls sharply, and stays 
small until it rises again near the phase boundary. Qual- 
itatively, this doping-dependence of the low-temperature 
specific heat coefficient agrees well with the ex peri mental 
results for high-T c superconducting materials [33 ■ 

V. CONCLUSIONS 

We have developed a position-space renormalization- 
group approximation to study the superfluid weight of 
the three-dimensional tJ model. Our results indicate 
that optimal hole doping for this system occurs in the 
density range of the r phase, (rii) w 0.63-0.68, where 
n s /m* reaches a local maximum. While the superfluid 
weight drops off sharply in the overdoped region, the den- 
sity of free carriers, proportional to the kinetic energy, 
remains approximately constant, as seen experimentally 
in high-T c materials. From calculations of the specific 
heat coefficient 7, we see clear evidence of a gap in the 
excitation spectrum for the r phase. Earlier renormaliza- 
tion group studies 0,0] had suspected that the t phase 
corresponds to the superconducting phase of high-T c ma- 
terials, and this idea was reinforced when an analogous 
phase was discovered in the Hubbard model Our 
present results justify this suspicion. 
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